Two Lines or Not Two Lines? That is the Question of Gamma Ray Spectra 
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Lines in the spectrum of cosmic gamma rays are considered one of the more robust signatures of 
dark matter annihilation. We consider such processes from an effective field theory vantage, and 
find that generically, two or more lines are expected, providing an interesting feature that can be 
exploited for searches and reveal details about the underlying theory of dark matter. Using the 
130 GeV feature recently reported in the Fermi-LAT data as an example, we analyze the energy 
spectrum in the multi-line context and find the data to be consistent with a single 77 line, a single 
"/Z line or both a 77 and a 7Z line. 

PACS numbers: 



INTRODUCTION 

Nonbaryonic dark matter is now a crucial element 
in the picture of the cosmology of the early Universe. 
And yet, its role in the framework of particle physics 
has remained elusive. Discovery of any kind of non- 
gravitational interactions of dark matter, through obser- 
vation of its production at high energy accelerators, its 
direct scattering with heavy nuclei, or its annihilation 
is an area of major experimental activity. Any of these 
observations would establish dark matter as an exotic 
particle, and would give insights into its nature. 

In the search for WIMP annihilation, gamma rays pro- 
vide a promising window. Gamma rays produced in the 
galaxy do not typically scatter on their way to the Earth, 
providing a handle from the morphology of their origin. If 
dark matter annihilates into quarks (or any particle with 
large decay branching ratios into quarks, such as W , Z, 
and Higgs bosons), the resulting spectrum of gamma rays 
tends to be rather soft, arising from the eventual decays 
of 7t s produced in the hadronic showers, and with a cut- 
off at the mass of the WIMP. These continuum signals 
are difficult to extract from the (often unknown) astro- 
physical backgrounds, and so to date searches have been 
most efficient when observing regions of the sky which are 
largely background free, such as dwarf spheroidal galax- 
ies in. 

If dark matter annihilates into charged particles, it 
must also be able to annihilate directly into two-body 
final states including a photon. Such processes are me- 
diated by loops, and thus are suppressed compared to 
the continuum annihilation. Their power comes from the 
feature that they produce a photon with a well-defined 
energy given (for the process XX ~^ 7^0 by 



E~ = m x 
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where My is the mass of the second annihilation prod- 
uct. For the case where Y = 7, the line occurs at an 
energy equal to the mass of the WIMP itself, E 1 = m x . 



Given this striking feature, the search for gamma ray 
lines has become a standard item on the menu of searches 
for WIMP annihilation using Fermi-LAT data [10. 

While it is possible for instrumental effects or more 
prosaic astrophysical processes 0] to mimic a bump in 
the gamma ray spectrum, a line remains one of the most 
compelling prospects for the indirect detection of dark 
matter annihilation. Motivated by the recent tentative 
indication that there may be such a feature at an energy 
around 130 GeV 0,0] with a relatively large (rough 1/10 
of the thermal expectation) cross section, and consistent 
with originating close to the galactic center 0-0 we ex- 
plore the generic properties that one might expect in a 
theory which can produce strong line signals. We use 
the Fermi-LAT data as analyzed in Ref. [6| to illustrate 
how one may dissect a putative line signal both to lend 
strength to its origin from dark matter annihilations, as 
well as to learn something about the details of the theory 
of dark matter. 



THE THEORY SPACE OF 7 RAY LINES 

Boiled down to its essence, the process XX ~^ lY re- 
sults from an amplitude involving a loop of charged par- 
ticles which also couple to Y. The charged particles in 
the loop could be either exotic states, or part of the Stan- 
dard Model, or (as is typical) a mixture of the two. For 
the current discussion, rather than wed ourselves to any 
specific UV-complete theory, we work in an effective the- 
ory framework and discuss operators in the effective ac- 
tion allowing WIMPs to annihilate into two particle final 
states, one of which is a 7-ray. 

The operator description is only appropriate to de- 
scribe theories for which the momentum transfer is 
smaller than the masses of any of the mediators which 
have been integrated out. For annihilation, the momen- 
tum transfer is ~ m x , so this restriction boils down to 
the requirement that all of the charged loop particles are 
much heavier than the WIMP itself (but one can enlarge 
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the effective theory to capture cases where some of the 
loop particles are heavier than the WIMP, and some are 
SM particles, see Ref. [§] for an example and Ref. @ for 
some related discussion). 

In constructing the operators in the EFT language, we 
consider both scalar and fermionic WIMPs. We iden- 
tify the leading operators of each type, and ignore higher 
mass dimension operators which are presumably further 
suppressed by the heavy mediator masses. We work in a 
description where the SU{2)l xU(1)y gauge symmetry is 
realized manifestly. In counting the dimension of opera- 
tors, this choice is equivalent to the assumption that the 
charged particles which have been integrated out have 
masses largely independent of electroweak symmetry- 
breaking. 

The natural ingredients from which operators are built 
are the field strengths of the hypercharge and SU(2) 
gauge fields, and W^ vl the Higgs doublet $, and 
the dark matter field \ (which we take to be either a 
scalar or spin- 1/2 fermion). We will build operators up 
to dimension-7 out of these ingredients, focusing on op- 
erators which produce at least one photon. We clas- 
sify each operator according to what type(s) of line pro- 
cesses) it mediates, including 77, 7Z [1CH13| . and/or 7/1 
[3, [H| (where h is the SM Higgs boson) 1 . We further 
distinguish operators leading to velocity-suppressed or - 
unsuppressed rates into gamma ray lines depending on 
whether the leading term in the expansion of (av) for 
small relative WIMP velocity v is a constant or is pro- 
portional to v 2 . Since v ~ 10~ 3 in the Milky way halo, 
velocity-suppressed operators need much larger couplings 
in order to produce a visible signal to compensate their 
v 2 suppression. 

At dimension-4, for either a scalar or fermionic WIMP, 
the unique choice leading to coupling to the pho- 
ton requires that the dark matter particles themselves 
are charged under t/(l)y and/or SU(2)l leading after 
EWSB to couplings to 7 and Z. The cosmological bounds 
on such "milli-charged" WIMPs are very strong [l7j . 
leaving this possibility unlikely to produce a line feature 
that could be observable by any near future experiment. 

At dimension-5, there are two operators for a Dirac 
fermion 

XTTxB^ and xfxB^ (2) 

where = 1/4[7 M , 7"]. These operators correspond to 
a weak magnetic (electric) dipole moment for the WIMP 
and leads to unsuppressed annihilation into both 77 and 
7Z. 



1 We necessarily will miss annihilation into a photon together with 
a dark sector particle [T3 | , which would require extension of the 
effective theory to contain the second decay product. 



At dimension-6, there is a family of operators built out 
of the set of dimension-4 factors, 

[B^B^, W^W*"", B^B» u , W^W a ^} (3) 

multiplied by y 2 (or |yj 2 if \ is a complex scalar). These 
operators also lead to unsuppressed annihilation into 
both 77 and jZ. 

Similarly, at dimension-7, there is another family of 
operators built out of the same set §Ji§ of dimension-4 
factors, multiplied by 

XX or xi5X ■ (4) 

for either a Majorana or Dirac fermion. These opera- 
tors generate both 77 and "fZ annihilations which are 
unsuppressed (suppressed) for X75X (xx)- 

Finally, for a Dirac fermion, we can also have the 
dimension-7 terms involving tensor operators formed 
from xi^X multiplied by a factor from either the set, 

{b^B™, W a m W a ™}, (5) 

or 

{B^\ 2 , B^ 2 , $tW^T a $, tfW^T**}, (6) 

the set contained in ([5]) leads again to both 77 and jZ 
lines (unsuppressed), whereas the set contained in ([6]) 
leads to a single unsuppressed 7ft- line. 

At dimension-8, there is a very large number of opera- 
tors, which we will not catalogue exhaustively. We note, 
however, that at dimension-8, there are operators built 
out of the vector currents, J M = {xduX* ~ X*d^X) (f° r a 
complex scalar WIMP), 5 M = xi^X f° r a Dirac fermion 
WIMP, and = X7/i75X f° r either a Dirac or Majorana 
fermion WIMP. Any of them may be combined with a 
factor from the set, 

{b^&D^, B m &D a $, $tw; Q ™ Q $, 

&W£ a T a D a $} . (7) 

All operators in this set lead to both 7Z and 7/1 lines. 
Those constructed from J p and are w-suppressed, 
whereas is unsuppressed. These operators are nat- 
urally generated in models where the dark matter dom- 
inant communication with the SM particles is through 
exchange of a Z' boson (for examples of models where 
this is the case, see [Il^ lTcl - Eoj ]). The absence of a 77 line 
is naturally explained by the Landau- Yang theorem [2l[ , 
which forbids a neutral vector state from decaying into 
two photons. 

The simple exercise of cataloguing operators in the ef- 
fective field theory already reveals very interesting fea- 
tures. Every operator considered leads to two lines (77 
and jZ) or jZ and 7/1), with a simple prediction for 
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their relative intensities. The sole exception is the set of 
operators contained in (|6]) which exist only for a Dirac 
WIMP and lead to a single 7/1 line. 

In terms of the underlying picture in which line pro- 
cesses are mediated by charged particles running in a 
loop, the fact that there are multiple lines corresponds 
to the fact that such particles must be charged under 
SU(2)l and/or U(l)y and thus must couple to both the 
photon and the Z boson. Similarly, coupling to a Higgs 
boson often is accompanied by coupling to a longitudinal 
Z boson (but not always - if the coupling is to it 
implies interactions with one or two Higgs bosons, but 
always two longitudinal Z bosons). Presumably any re- 
alistic UV complete theory will generate more than one 
operator, allowing for the possibility of interference. In- 
terference will adjust the relative sizes of the two or three 
lines, but in the absence of fine tuning will not cancel one 
of them completely. 

Our results suggest interesting variations in the ex- 
perimental analyses of gamma ray lines. First, one may 
search for two (or three) lines whose energies are con- 
sistent with Eq. (JT|) for a single WIMP mass. Such an 
observation would be highly suggestive of dark matter 
annihilation, and less likely to be produced by astro- 
physical or instrumental effects. Of course, the ability 
to resolve two lines is very challenging for WIMP masses 
larger than ~ 150 GeV, because of the finite energy res- 
olution of the detector (AE/E ~ 10% at E ~ 100 GeV 
for the Fermi-LAT). Second, if there is a concrete obser- 
vation of a single line, will be highly suggestive of a Dirac 
fermion WIMP annihilating through one of a definite set 
of effective interactions, providing clues to the nature of 
the UV theory. 

THE FERMI 130 GEV FEATURE AS A CASE 
STUDY 

In order to explore how multiple lines could mani- 
fest themselves in realistic data, we analyze the observed 
7-ray spectrum of Fermi-LAT, including the feature at 
E~/ — 130 GeV 0, in the context of the multi-line the- 
ory. While it is premature to interpret this feature as 
dark matter annihilation, it provides an interesting case 
study which could even turn out to ultimately tell us 
something about dark matter. We use the regions of in- 
terest found to have largest significance, Reg3 and Reg4 
[!, Q and the ULTRACLEAN photon selection. While 
not presented here, we have also performed the analysis 
for the looser photon selection (SOURCE class), which 
yields very similar results. We will focus our analysis on 
the 77 plus jZ case; the extrapolation for 7Z plus 7/1 is 
straightforward. 

We follow the standard Fermi analysis procedure, eval- 
uating the relative likelihood of the background-only hy- 
pothesis (null) and the background-plus-signal (best) hy- 



pothesis using the test statistic (TS): 

TS = -21n^i (8) 

but with a likelihood C which includes both a power-law 
background model as well as terms for potential 77 and 
r)Z lines: 

£(£ 7 |V 77 ,V 7Z ,/3,a) = 

P (J) + V 77 / DM (£ 7 K) 

+ ^/dm(^K^1-0)) (9) 

The function /dm Inline) is a normalized double 
Gaussian function fit to the expected line shape for 
-Eiinc = 130 GeV as provided in Ref @. For other val- 
ues of the expected peak location i?n ne , the values of 
the Gaussian widths and means are treated as linearly 
dependent on the position of the expected peak. The pa- 
rameters iV 77 and N^z control the total yield from the 77 
and 7Z processes, respectively. The two terms describe 
the correlated 77 and 7Z contributions. 

For £best, the parameters {N 11 ,N 1 z,P,Oi) are floated 
to find the maximum likelihood value. For £ n uii, the 
yields (iV 77 , N 7 z) are fixed to zero and the background 
model parameters (j3, a) are floated to their best fit val- 
ues. The local statistical significance may be interpreted 
in the asymptotic regime [22j |. as a = VTS. 

Example fits are shown in Fig [TJ If the WIMP is as- 
sumed to have mass of 145 GeV, a "fZ process would pro- 
duce a line at E 1 — 130 GeV. In this case, the feature at 
E-y = 130 GeV can be interpretted as pure "fZ; any con- 
tribution from 77 would appear at larger E 7 , where no 
such feature appears. If, on the other hand, the WIMP 
is assumed to have a mass of 130 GeV, then a 77 process 
would produce a line at E~ = 130 GeV, explaining the 
feature 2 . In this case, however, there is room at lower 
values of _E 7 for contributions from a "fZ process. The 
two fits have approximately equal signifiance. 

Figure [2] shows a scan of WIMP masses and total yields 
(iV 77 + Njz), revealing the two regions of maximal sig- 
nificance, near m x — 130 and 145 GeV. At each point, 
the contribution from 77 is also shown. If mwiMP = 130 
GeV, then the interpretation is consistent with a large 
77 contribution, but cannot rule out some contribution 
from 7Z. If, however, m x = 145 GeV, then a pure 7Z 
interpretation is preferred. Figure [3] shows the same scan 
for Reg4, which has largely the same features. 



2 In Ref. it is claimed that in order to produce a peak at E~, — 
130 GeV after electromagnetic showering of the parent photon, 
a WIMP mass of ~ 145 GeV is required. We disagree with this 
statement. 
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FIG. 1: Sample spectrum predicted by dark matter annihila- 
tion into a single 7Z line (top panel) or two lines produced 
by 77 and 'yZ (bottom panel). Data from Ref. G\ is overlaid. 



FIG. 2: Top, local statistical significance of the signal ver- 
sus WIMP mass and total yield due to 77 and 7Z lines in 
the Reg3 region of interest defined in Ref. ||. Bottom, the 
contribution from 77 at each point in the plane. 



In Figures [4] and [5] the scans are performed as a func- 
tion of WIMP mass and the ratio of to 77 yields 
for Reg3 and Reg4. Note that while the region of max- 
imum significance is near m x — 130 GeV, it prefers 
N 1 z/N 11 > 0. Also shown in each figure is the value 
of iV 77 corresponding to the best fit for each point in 
(m x ,N lZ /N 17 ). 

The two regions show consistent features. The maxi- 
mum significance is consistent with either a pure 7Z or 
pure 77 scenario; the interpretation of the E 1 = 130 GeV 
line in the 77 scenario also allows for 7Z contributions 
at lower E y . In fact, the best fits have a non-zero frac- 
tion of 7Z (less than one) , but this preference is not very 
significant. Clearly more data would be very helpful in 
terms of sharpening this analysis in order to draw more 



firm conclusions from it. 

From Ref. @, the best fit value of the cross section 
(assuming a 77 interpretation) is about 10 -27 cm 3 /s 
~ 10~ 4 TeV~ 2 for dark matter distributed according to 
an NFW profile. From here one could compare with de- 
tailed calculations based on the operators in the effective 
field theory catalogue to determine a consistent param- 
eter space, but we leave such detailed comparisons for 
future work, and instead interpret such a target cross 
section schematically. 

Focusing as an example on any one of the dimension-6 
operators for scalar WIMPs, we normalize the operator 
as aa x /M 2 , where a = e 2 /(47r) is the electromagnetic 
coupling and a x = g 2 /(47r) represent (unknown) cou- 
plings in the dark sector. This choice of normalization is 
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FIG. 3: Top, local statistical significance of the signal ver- 
sus WIMP mass and total yield due to 77 and 7Z lines in 
the Reg4 region of interest defined in Ref. [g|. Bottom, the 
contribution from 77 at each point in the plane. 

consistent with the operator being generated at one loop, 
with M playing the role of the mass(es) of the particles 
in the loop. Obviously, this implementation is subject to 
unknown numerical factors such as the number of species 
contributing inside the loop, as well as factors associated 
with their spins, etc. The idea is to get a very rough 
sense for the mass scale M of the loop particles, given 
the target cross section of the Fermi feature. 

Our simple estimate indicates that provided there is 
no velocity-suppression, 

M ~ i/a^ 150 GeV. (10) 

This is an interesting result. If the dark sector is strongly 
coupled (a x ~ 1), the loop particles should have masses 
in the range of 150 GeV, safely above the LEP bound 



Reg3 ULTRACLEAN 




FIG. 4: Top, local statistical significance of the signal versus 
WIMP mass and ratio of yields in the 77 and 7Z lines in 
the Reg3 region of interest defined in Ref. ||. Bottom, the 
contribution from 77 at each point in the plane. 



of about 100 GeV, but low enough that the LHC has an 
opportunity to observe them through electroweak pro- 
duction. For weaker a x , the mass must be lower to com- 
pensate, rapidly coming into conflict with the LEP bound 
for a x < 0.5. 

For a velocity-suppressed operator, the target mass is 
of the order M ~ y^hc 5 GeV, far enough below the 
LEP bound that not even a strongly coupled dark sector 
would be able to reconcile the two. However, it is worth 
mentioning a few provisos to this statement. For exam- 
ple, one way in which the EFT could spectacularly break 
down would be when there is an additional dark sec- 
tor state which appears in the s-channcl for annihilation. 
Very large enhancements are possible in this case, de- 
pending how close to on-shell the resonance is for WIMP 
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The effective theory illustrates a fascinating feature - 
the operators which give rise to one gamma ray line, typ- 
ically also give rise to two. For a scalar or Majorana 
WIMP, every operator considered produces either 77 and 
"fZ, or jZ andjh, and the intensities of each line are cor- 
related for a given operator. Multiple lines are a generic 
feature, and one that can be used to improve searches 
in data from gamma ray observatories, or help match to 
specific UV complete theories once a discovery is made. 
For a Dirac WIMP, one class of operators provides an 
exception to the multiple-line rule, producing a single 7/1 
line. Nonetheless, observation of a single line provides 
very specific information about the nature of the theory 
of dark matter. 

Using the recent observation of a feature at 130 GeV 
in the Fermi-LAT data, we analyze the data in a multi- 
line context, and find that there is a very mild prefer- 
ence for contribution from two lines, though uncertainties 
are large. Should this feature persist and not ultimately 
prove to be instrumental or astrophysical in nature, more 
data should help sharpen this analysis and make more 
concrete statements. 
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FIG. 5: Top, local statistical significance of the signal versus 
WIMP mass and ratio of yields in the 77 and 7Z lines in 
the Reg4 region of interest defined in Ref. [f|. Bottom, the 
contribution from 77 at each point in the plane. 

annihilation. In addition, if a large multiplicity of species 
contribute to the line annihilation, the amplitude will 
grow with the number. For a rather extreme multiplicity 
of ~ 500, a w 2 -suppressed annihilation would be consis- 
tent with the LEP bound for a x ~ 1. 



OUTLOOK 

Annihilation of dark matter into a two body final state 
containing a photon provides a striking signature, and is 
one of the most promising prospects for an indirect de- 
tection of dark matter. In this article, we have explored 
some generic features of gamma ray lines using an effec- 
tive theory framework. 
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